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The magnetic field induced phase separation 
in a model of a superconductor with local electron pairing 
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We have studied the extended Hubbard model with pair hopping in the atomic limit for arbi- 
trary electron density and chemical potential and focus on paramagnetic effects of the external 
magnetic field. The Hamiltonian considered consists of (i) the effective on-site interaction U and 
(ii) the intersite charge exchange interactions /, determining the hopping of electron pairs between 
nearest-neighbour sites. The phase diagrams and thermodynamic properties of this model have been 
determined within the variational approach (VA), which treats the on-site interaction term exactly 
and the intersite interactions within the mean-field approximation. Our investigation of the general 
case shows that the system can exhibit not only the homogeneous phases: superconducting (SS) 
and nonordered (NO), but also the phase separated states (PS: SS-NO). Depending on the values 
of interaction parameters, the PS state can occur in higher fields than the SS phase (field-induced 
PS). Some ground state results beyond the VA are also presented. 

PACS numbers: 

71.10.Fd — Lattice fermion models (Hubbard model, etc.), 
74.20.-z — Theories and models of superconducting state, 

64.75.Gh — Phase separation and segregation in model systems (hard spheres, Lennard-Jones, etc.), 
71.10.Hf — Non-Fermi-liquid ground states, electron phase diagrams and phase transitions in model systems 



I. INTRODUCTION 

Recently, there has been much interest in supercon- 
ductivity with very short coherence length. This inter- 
est is due to its possible relevance to high temperature 
superconductors (the cuprates, doped bismuthates, iron- 
based system, fuUerenes, etc.; for a review, see [T}{5] and 
references therein). Also the phase separation (PS) phe- 
nomenon involving superconducting (superfluid) states is 
very current topic, because it can play crucial role deter- 
mining behaviour in many real compounds |5HT5] and 
fermions on optical lattices 



The Hamiltonian considered has the following form: 
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where fii = Ea^J'r, fii^ = cj^h^, p] ^ {p^ y ^ c^^c-^. 
B = gpsHz is external magnetic field and 
— (1/2) (rii^ — fiii) is z-component of the total spin at 
i site, denotes the creation operator of an electron 
with spin (t=^,^ at the site i, which satisfy canon- 
ical anticommutation relations: {cia,cf^,} = SaS, 



is the Kronecker 



In present work we will study paramagnetic effects of 
magnetic field (Zeeman term) in a model which is a gen- 
eralization of the standard model of a local pair super- 
conductor with on-site pairing (i. e. the model of hard 
core bosons on a lattice [T] [5TJ |21]) to the case of fi- 
nite pair binding energy. Such analysis of paramagnetic 
effects is important, in particular, for unconventional su- 
perconductors, for which the temperature dependence of 
the upper critical field has positive curvature PSI 12^ and 
does not saturate even at genuinely low temperature |24| . 
Also the pseudogap is destroyed by sufficiently high field 
[25J. Moreover, recently the possibility of the magnetic 
field induced phase separation (SS/NO) has been found 
for the d = 2 and d = 3 dimensional spin-polarized at- 
tractive Hubbard model ^26j as well as for the continuum 
fermion model in d = 2 |27| . 
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{cia,Cja'} = {cZr,Cj^'} = 0, whcrc S, 
delta. j") indicates the sum over nearest-neighbour 
sites i and j independently. U is the on-site density 
interaction, / is the intersite charge exchange interaction 
between nearest neighbours, p is the chemical potential, 
connected with the concentration of electrons by the 
formula: n = {1/N) {f^i) i with < n < 2 and TV is 
the total number of lattice sites. 

In this paper, we treat the parameters U and / as 
the effective ones, assuming that they include all the 
possible contributions and renormalizations like those 
coming from the strong electron-phonon coupling or 
from the coupling between electrons and other elec- 
tronic subsystems in solid or chemical complexes jl]. 
In such a general case arbitrary values and signs of 
U and / are important to consider. Formally, / is 
one of the off-diagonal terms of the Coulomb interac- 
tion lij — — (l/2)(M|e^/r|jj) ^28j i29j, describing a part 
of the so-called bond-charge interaction, and the sign of 
the Coulomb- driven charge exchange is typically nega- 
tive (repulsive, / < 0). However, the effective attrac- 
tive interaction of this form (/ > 0) is also possible [301 - 
I32| . In particular, it can originate from the coupling 
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of electrons with intersite (intermolecular) vibrations via 
modulation of the hopping integral |30j , or from the on- 
site hybridization term in generalized periodic Anderson 
model PT1I52]. 

In the analysis of model ([T]) at T > we have adopted a 
variational approach (VA) which treats the on-site inter- 
action U exactly and the intersite interactions / within 
the mean-field approximation [31 HI US]- Moreover, at 
T = exact results for d = 1 and results obtained within 
random phase approximation (RPA) (for d = 2 and d = 3 
lattices) are presented. 

Model ([ij has been investigated intensively in the ab- 
sence of the external magnetic field only. The first anal- 
ysis of the phase diagram of model (jlj have been per- 
formed by Bari |M| and Ho and Barry [3S] using the 
variational method in order to examine the instability of 
the Mott insulator to superconductivity mostly for the 
special case of the half-filled band {n — 1). The effects of 
diagonal disorder on the critical temperature for U — 
and n = 1 have been also determined |36j . arriving at 
a satisfactory qualitative interpretation of quite a num- 
ber of different experiments in amorphous superconduc- 
tors. Within the VA the phase diagrams of model ([T]) 
as a function of the electron concentration n for i? = 
have been investigated in (3] 0] HI] • The stability condi- 
tions of states with phase separation for B = have been 
discussed in [3] only. 

In this paper, we investigate model Hamiltonian ([T]) 
for arbitrary /i and arbitrary n at T = and finite tem- 
peratures. We focus on the effects of external magnetic 
field in the system. Our investigation of the general case 
finds that, depending on the values of the interaction pa- 
rameters and the electron concentration, the system can 
exhibit homogeneous SS and NO phases as well as the PS 
between them. Transitions between various states and 
phases can be continuous and discontinuous, what im- 
plies existence of tricritical points on the phase diagrams. 
We present detailed results concerning the evolution of 
the diagrams as a function of external field B, interaction 
parameters, p, and n and discuss representative thermo- 
dynamic properties of the system. The results obtained 
in this work can be useful for the description of systems 
with local pairing in the magnetic field. They are also im- 
portant as a test and a starting point for a perturbation 
expansion in powers of the hopping tij and as a bench- 
mark for various approximate approaches analyzing the 
corresponding finite bandwidth models. 

The paper is organized as follows. In section [IT] we de- 
scribe the VA method. Section [iHjisde voted to the study 
of the phase diagrams: sectioii~ |lII A| includes results at 
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represen- 



T = (VA, exact, RPA) whereas in section 
results at T > are presented. In section 
tative thermodynamic characteristics are evaluated and 
discussed. Finally, section [V] contains conclusions and 
supplementary discussion. 



II. THE VARIATIONAL METHOD 

Within the VA the on-site interaction term is treated 
exactly and the intersite interactions are decoupled 
within the mean-field approximation (site-dependent): 

PtpJ ^ (Pt) Pj + (Pj) Pt - {pt) (Pj) ■ (2) 
A variational Hamiltonian has the following form: 

Ho 
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+ X*A^ + X^K~Bst], (3) 

where Xi = Y.jJ.ihj^j, A* = {pf) and = (n,). 
can be diagonalized easily and a general expression for 
the grand potential J7 in the grand canonical ensemble in 
the VA is 



-^In^Tr 



{Tr [exp(-/3i/o)] } 



where /3 = l/{kBT). The average value of operator A is 
defined as 
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Tr_B means a trace of operator B calculated in the Fock 
space. 

Assuming no spatial variations of the order parameter 
the grand potential per site obtained in the VA for model 
([ij is given by: 

uj{p) = Vl/N = + 2/o|A|2 - ln(2Z), (4) 

where 

Z = cosh f/3v/Ai2 +4|/oA|2) + exp {pU/2) cosh (/3B/2), 



p = pi- C//2, /o = zl, and A* = {l/N) z de- 

notes the number of nearest neighbours. The free energy 
per site f — uj + j.m is derived as 

f{n) = fl{n - 1) + (U/2)n + 2/o|A|2 - f]-^ ln(2Z). (5) 

In the absence of the field conjugated with the 
SS order parameter (A) there is a symmetry be- 
tween / > (s-pairing) and / < (?7-pairing, r/S, 
A^5 — J2i 6xp (iQ ■ Ri){pJ), Q is half of the smallest 
reciprocal lattice vector) cases for model Q, which ne- 
glects single particle hopping {tij = 0), thus we restrict 
ourselves to the / > case only. 

The condition for electron concentration and a mini- 
mization of w (or /) with respect to the superconducting 
order parameter |A| lead to the following self-consistent 
equations (for homogeneous phases): 



/Isinh (^^/i2+4|JoA|2 



ZVm2 + 4|/oA|2 
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sinh (^/3Va*'+4|/oA|2 
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FIG. 1. Phase diagrams at T = 0: (a) as a function of /i//o, and (b) as a function of n. Solid, dotted, dashed, and dashed-dotted 
Hnes denote the boundaries derived within VA (exact for d — >■ +c»), d = 1 (rigorous results), d = 2 (RPA, SQ lattice), and d = 3 
(RPA, SC lattice), respectively. At half-filling (n = 1) the NO (Mott state) is stable above the end of PS-SS lines. The NO-SS 
transition for p./Io = —1 (panel (a)) is second order, all other transitions between homogeneous phases are first order. The 
transitions to the PS state are "third order". 



Equations ([6])-(|7]) are solved numerically for T > and 
we obtain |zXj and n when fj, is fixed or | A| and n when n is 
fixed. The superconducting phase (SS) is characterized 
by non-zero value of |A|, whereas |A| = in the non- 
ordered (normal) phase (NO). 

The magnetization of homogeneous phases can be sim- 
ply derived as 



dn 

dB 



2Z 



exp(/3t//2) sinli(/3B /2) . (8) 



It implies that the magnetization in both SS and NO 
phases is nonzero for any B ^ and T > 0. In the ho- 
mogeneous phases the double occupancy per site, defined 
as D — J2i {fii^nii), has within VA the following form: 



D = 



l-^exp(/3(7/2) 
nZ 



(9) 



We also introduce the concentration of locally paired elec- 
trons Tip — 2D and the ratio Up/n = 2D/n. Notice that 
D is different from the condensate density (a fraction of 
pairs in the condensate) no = |(/5^)p. 

Phase separation is a state in which two domains with 
different electron concentration: n+ and n_ exist in the 
system (coexistence of two homogeneous phases). The 
free energies of the PS states are calculated in a standard 
way, using Maxwell's construction (e. g. [11|37]). In the 
model considered only one type of PS states can occur, 
which is a coexistence of SS and NO phases. 

In the paper we have used the following conven- 
tion. A second (first) order transition is a transition 
between homogeneous phases with a (dis-)continuous 
change of the order parameter at the transition temper- 
ature. A transition between homogeneous phase and the 
PS state as a function of n is symbolically named as 
a "third order" transition |4l [37l |55]. This denotation 
should not be misled with the Ehrenfest's notation of or- 
der of transitions between homogeneous phases. At this 



transition a size of one domain in the PS state decreases 
continuously to zero at the transition temperature. One 
should notice that the order parameter for "third order" 
transitions is the concentration difference n_|_ — n_ (not 
A, which is the order parameter in one domain) and its 
change is discontinuous at transition temperature. Such 
transitions are present if the system is considered for fixed 
n and they are associated with first order transitions at 
fixed /2. 

All phase transition boundaries, necessary to construct 
the complete phase diagram within VA, have been ob- 
tained numerically by self-consistent solving of ([6])-([7]) 
and comparing grand potentials uj of homogeneous phases 
(if /2 is fixed), or free energies / - including energies of 
PS states - if n is fixed. 



III. PHASE DIAGRAMS 

The diagrams obtained are symmetric with respect to 
half-filling (n = 1) because of the particle-hole symmetry 
of Hamiltonian ([ij , so they will be presented only in the 
range /2 < and < n < 1. 



A. The ground state 

The ground state energy of the SS phase within the 
VA is derived as fss = (l/2)C/n - (l/2)|/o|n(2 - n) with 
|Ap = (l/4)n(2 - n) and m = 0. For the NO phase at 
T = 0, if n is fixed, one obtains fjyo = —{l/2)Bn (with 
2m = n for i? ^ and 7i < 1). If /i is fixed, one has (i) 
^No{P^) = (n = 0) and (ii) uj^oiP) ~ —fi — U/2 — B/2 
(2rn = n = l, NO - Mott). Notice that for the NO 
phase the VA gives the rigorous results for the / 
and u!. The energy of the PS state is obtained by 
Maxwell's construction. In the NO domain njvo — 1 
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FIG. 2. Finite temperature phase diagrams for U/Io — (a), U/Io = 0.75 (b), and U/Io = 1.25 (c) plotted as a function of 
fi/Io and B/Iq. Solid and dotted lines indicate second order and first order transitions, respectively. The tricritical point line 
is denoted by dashed-dotted line. 




n n n 



FIG. 3. ksT/Io vs. n phase diagrams for U/Iq — 0, B/Iq = 0.99 (a) (corresponds to form (ii)), U/Iq = 0.5, B/Io = 0.5 (b), and 
U/Io ~ 1, B/Io = 0.25 (c) (corresponds to form (iii)). Details in text. Solid and dashed lines indicate second order boundaries 
and (named symbolically) "third order" boundaries (i. e. between the PS state and homogeneous phases), respectively. T 
denotes the tricritical point. 



(Mott phase with m = 1/2 for i? 7^ 0) , whereas concen- 
tration in the SS domain is nss = 1 ^ \/ {U + B) / Iq — 1 
{1 < {U + B)/Io < 2). It corresponds to the first order 
SS-NO boundary on the U/Iq vs. p,/Io diagram deter- 
mined by equation (p,/ 1^)'^ + 1 = U / Iq. In the NO phase 
m = for B = 0. 

It is also possible to obtain some results at T = be- 
yond the VA (being rigorous for d = 00) in dimensions 
d = 1, 2, 3 by the decomposition of the eigenspace of H 
into sectors specified by the parity of the occupation num- 
ber at each site, similarly as it has been done in [3] for 
B — 0. For d = 1, by making use of the exact results 
for the ground state of the d = 1 XY model in transverse 
field one can obtain the energy for the SS phase at T = 0. 
In dimensions 1 < d < +00 the self-consistent random 
phase approximation (RPA) is a reliable approach. It has 
been proven to be a very good approximation scheme in 
problems of quantum magnetism and it fully takes into 
account quantum fiuctuations, which can be of crucial 
importance for the considered system for d < 3. The 
resulting diagrams are shown in figures [T] For fixed 
n, the first order SS-NO transition line in the P/Iq vs. 
{U + B)/Io plane is replaced by the PS region bounded 
by two critical {U + B)/Iq values, the lower decreasing 
with n and the higher independent of n. In finite di- 
mensions due to quantum fiuctuations the regions of the 



homogeneous SS phase occurrence are extended in com- 
parison with the VA results. Moreover, at T = in d = 1 
only a short-range order occurs in contrary to the VA and 
RPA in d > 2, where the long-range order is present. 



B. The finite temperatures (within the VA) 

Let us discuss now the finite temperature phase 
diagrams as a function of ft (figure [2|. For 
U/Io < (2/3) ln(2) the tricritical point T, connected with 
a change of the transition order, appears only for B > 0. 
For any \p.\/Io < 1 at temperatures above T-point the 
SS-NO transition is of second order (figure |2^). For 
(2/3) ln(2) < U /Iq the line of tricritical points starts from 
the B = 0-plane. If (2/3) ln(2) < U/Iq < 1 the SS phase 
with n = 1 can occur at finite temperatures (figure ^p) , 
whereas for 1 < U/Iq < 2 the SS phase can appear on 
the diagram only for p- ^ (figure In the range 

2 < {U + B)/Io < +00 only the NO phase is stable at 
any B > and T > 0. For half-filling the T-point (if it 
exists) is located at ksT/Io = 1/3 (cf. also figure [5]). 

A few examples oi ksT / Iq ys. n phase diagrams are 
presented in figure [3l The dia grams, depending on given 
values of U/Iq and B/Iq, can be one of the following 
three forms only: 
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FIG. 4. B/Iq vs. kBT/Io phase diagrams for U/Iq = 1.25 and n = 0.1, 0.25, 0.75 (corresponding to forms (A)-(C) of the phase 
diagrams, respectively, defined in text). Denotations as in figure [s] 



TABLE I. Occurrence of the phase diagrams forms (/o — 1). 
These forms are defined in text referring to figures |3] and |4] 



Form 


(7 = 


U = 0.75 


U = 1.25 


(i) 


< B < 0.88 






(ii) 


0.88 < B < 1 


< B < 0.25 




(ui) 


1< B < 2 


0.25 < B < 1.25 


< B < 0.75 


(A) 


< n < 1 


< n< 0.3 


0<n< 0.11 


(B) 




0.3 < n < 1 


0.11 < n < 0.5 


(C) 






0.5 < n < 1 



(i) only homogeneous phases occur on diagrams at any 
T and the SS-NO transition is of second order, the 
transition temperatures increase with increasing n, 

(ii) the PS state appears only at T > 0, e. g. figure |3|i 
(the SS-NO first order fine ends at T > on the 
ksT/lQ vs. fi/Io diagram for the same model pa- 
rameters [///q and B/Iq), 

(iii) the PS region extends from T = (for 
1<{U + B)/Ia < 2), e. g. figure ^ (the SS-NO 
first order line ends at T = on the ksT/If) vs. 
fi/Io diagram). 

Figure pp shows the limiting case between form (ii) and 
form (iu). Ranges of B/I^ in which forms (i)-(iii) occur 
for three particular values of U/Iq are collected in the 
upper part of Table |T1 

One can distinguish three different forms of B/Iq vs. 
ksT/Io phase diagrams obtained for a given fixed n 1 
(and U/Iq), which are shown in figure H (numerical val- 
ues for U/Io = 1.25). For U/Iq < (2/3)ln(2) only the 
form (A) (shown in figure J4a) occurs for any < n < 1 
(cf. bottom part of Table In such a case the tricriti- 
cal point T exist at T > 0. At higher temperatures and 
lower fields the SS-NO transition is continuous (cf. also 
figure [2|3). At lower temperatures and higher fields the 
SS and NO phases are separated by the PS state. In the 
range (2/3) ln(2) < J7//o < 1 form (A) is realized for suf- 
ficiently small n, whereas for higher n form (B) appears, 
where the regions of the SS and NO phases occurrence on 
the diagrams are separated by the PS state for any field 




0.0 0.2 0.4 0.6 0.8 

B/r 



FIG. 5. feflT/Jo vs. B/Iq phase diagrams for n = 1 and vari- 
ous U/Io (as labelled). Denotations as in figure [5] 



and temperature (figure |4p) . For 1 < U/Iq < '2 and n 
sufficiently close to half-filling only the PS state appears 
on the diagrams (absence of the homogeneous SS phase, 
form (C) of phase diagrams, figure |4|;). Ranges of n in 
which forms (A)-(C) occur for three particular values of 
U/Iq are given in the bottom part of Table |l] 

The resulting phase diagrams for n — I and various 
U/Io are shown in figure [s] In such a case only homoge- 
neous phases exist on the diagram. The SS-NO transition 
with increasing T is of second order for ksT/Io > 1/3 
(it occurs only if U/Iq < 0.46) and of first order for 
keT/lQ < 1/3 (and any U/Iq < 1). 

Concluding this section, the possible sequences of tran- 
sitions with increasing temperatures and the transition 
orders of them for the system at fixed n are listed below: 

(i) SS— ?>N0: second order for n ^ 1 and second order 
or first order for n = 1, 



(ii) PS— >N0: "third order", it can take place only for 

(iii) SS-^PS^NO: both "third order", it can take place 
only for n ^ 1. 
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FIG. 6. Temperature dependencies of (a) the superconducting order parameter |A| (solid lines) and the magnetization 3m 
(dashed lines) (b) the ratio Up/n and (c) the specific heat c/fes plotted for: p./Io = —0.1, U/Io ~ —0.5 and B/Io = 0.5 (A), 
fl/Io = -0.1, U/Io = 0.75 and B/Io = 0.1 (B) as well as for: fi/Io = -0.8, U/Io = 0.75 and B/Io = 0.1 (C), fi/Io = -0.8 
U/Io = 1.25 and B/Io = 0.1 (D). 



IV. THERMODYNAMIC PROPERTIES 

In this section we present several representative depen- 
dencies of the thermodynamic characteristics for fixed 
model parameters. In particular, for fixed p, (figure [6|, 
one can single out two limiting types of thermodynamic 
behaviour near transition temperature Tss' (i) the lo- 
cal pair regime and (ii) the pair breaking regime. In be- 
tween, there is a crossover between the two regimes. Let 
us stress that in our model the single particles do coexist 
with pairs at finite temperatures (except U — )• — oo) and 
at T > the concentration of paired electrons Hp is al- 
ways smaller than n and there exist finite concentration 
of single particles: n — rip in the system. It does modify 
the phase transitions, which properties in the local pair 
regime (np{Tss) < n) are different then those in the pair 
breaking regime, where np{Tss) is substantially smaller 
than n (cf. figure [gJ^). 

For large on-site attraction the concentration of locally 
paired electrons Up = 2D exhibits no sharp feature as the 
temperature is lowered through Tss (line A in figure |6|3). 
The number of non-paired electrons at Tss is negligible 
and the transition is to the state of dynamically disor- 
dered pairs (only for \U\/Io 3> 1 and U < 0, the local 
pair regime). In the second limit, for on-site repulsion 
U < 2/o, Hp has a sharp break at Tss E^nd a substan- 
tial fraction of single particles can exists both below and 
above Tss- We call this the pair breaking regime. As 
temperature is lowered, the condensate grows both from 
a condensation of pre-existing pairs and from binding and 
condensation of single particles. For small binding ener- 
gies, if n ^ 1 {fl/Io « —1), there will be essentially no 
pre- formed pairs at Tss (lines C and D). 

Obviously a more realistic description of coexistence of 
itinerant electrons and local pairs in particular real ma- 
terials would be obtained using multicomponent mod- 
els, especially mixed boson-fermion model (see e. g. 

Notice that non-zero value of Up does not imply that 
local pairs are in coherent state and even significant val- 
ues of rip/n are possible in the NO phase. In the limit 



T +00 Up increases to Up/n — 0.5 (each of four states 
at a given site can be occupied with equal probabil- 
ity). The condensate density (which can be approxi- 
mated as no « |Ap at least for n <ti 1, uq ^ D) van- 
ishes for T>Tss, but the doubly occupied states are 
still present above Tss {D 0). 

In figure|6^ the temperature dependencies of the super- 
conducting order parameter A are presented (solid lines) , 
where one can see clearly the discontinuous change of the 
order parameter (lines B and D). The other lines corre- 
spond to second order transitions. The dependencies of 
magnetization m are denoted by dashed lines (lines for 
cases B and C are not distinguishable in the SS phase). 
One can notice that m in SS phase is strongly reduced 
(to = at T = 0) and increase with increasing T, whereas 
in the NO phase m decreases with increasing T and 
TO = for T — +00 (for —2/2 > U + B magnetization 
TO can increase with T near above Tssi e. g. cases C and 
D). 

Finally, let us briefiy summarize the behaviour of the 
specific heat at constant volume c — — T ^fj^^^ (fig- 
ure |6j;) . The NO phase is characterized by the relatively 
broad maximum in c connected with continuous changes 
in a short-range electronic ordering (in higher tempera- 
tures, not shown in figure |6|. The narrow peak in c{T) 
is associated with the first order transition, while the 
A-shape behaviour is typical for the second order transi- 
tion. 

The behaviour of thermodynamic parameters in PS 
states for the model with B ^ is similar to that for 
B = 0, which was widely discussed in section 5 of [4. 

V. CONCLUSIONS AND SUPPLEMENTARY 
DISCUSSION 

In this paper, we have studied the paramagnetic effects 
of external magnetic field on a simple model of a super- 
conductor with very short coherence length (i. e. with the 
pair size being of the order of the radius of an effective 
lattice site) and considered the situation where the single 
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particle mobility is much smaller than the pair mobility 
and can be neglected. In the model considered the pair 
binding energy is given hy Ei, = —U + 2Iq, whereas the 
pair mobility ^ / (for B — 0) and the critical magnetic 
field is proportional to the pair binding energy Eb (at 
least in the low concentration limit) [3_. 

Let us summarize important conclusions of our work. 

(i) For fixed /i, if for B = the SS-NO transition is of 
second order, the increasing magnetic field changes 
first the nature of phase transition from a continu- 
ous to a discontinuous type and next it suppresses 
superconductivity. 

(ii) In definite ranges of n and temperature the mag- 
netic field stabilizes the phase separation state 
SS-NO (field-induced PS). 

(iii) For fixed n, one can distinguish four different struc- 
tures of phase diagrams at T > 0, as illustrated in 
figure [4] (n 7^ 1) and figure [5](n = 1). 

These behaviours are associated with the presence of the 
tricritical point on the phase diagrams. 

For U —00 (states with single occupancy are ex- 
cluded and only local pairs can exists in the system) 
the model is equivalent with the hard-core charged boson 
model on the lattice (3] dll 112 and in such a case there 
is no paramagnetic effects oi B. In the local pair limit 
(for sufficiently strong on-site attraction and small B) 
the charge exchange / and on-site U < cooperate and 
the second order SS-NO transition is associated with the 
transition to a state of dynamically disordered local pairs. 
The opposite regime, i. e. the pair breaking limit, is re- 
alized for substantial values U and B {{U + B)/Iq ^ 1). 
In this limit the transition is determined by pair break- 
ing excitations (driven by both U and B, which destroy 
the electron pairs) and there are essentially no preformed 
pairs close to transition temperature (T > Tc). In gen- 
eral, for fixed U and p, (or n) the regions of ordered states 
occurrence are reduced by increasing B. 

Within the VA the on-site U term is treated exactly. 
Thus, the major conclusions of our paper concerning the 
evolution of the properties of the system with U are re- 
liable for arbitrary U. The derived VA results are exact 
in the limit of infinite dimensions d — > -l-oo, where the 
mean-field approximation treatment of the intersite in- 
teraction / term becomes the rigorous one. Moreover, 
the VA yields exact results (in the thermodynamic limit) 
for lij of infinite range (/,y — {1/N)I for any {i,j)) re- 
gardless of the dimensionality (d < 00) of the system f35]. 
However, for short range lij in finite dimensions due to 
quantum fluctuations the regions of the homogeneous 
SS phase occurrence are extended in comparison with 
the VA results (cf. flgure[T|). Also, in ID-chain only a 
short-range order occur at any T > 0, whereas in d = 2 
a long-range order can exist at T = 0, while at T > the 
Kosterlitz-Thouless transition is only possible. 

The presence of the hopping term j ^tijcf^Cja 
breaks a symmetry between the / > (favoring SS) and 



/ < (favoring rjS) cases. For tij 7^ model ([ij is called 
the Penson-Kolb-Hubbard model |3[0j. In general for 
tij ^ the phase diagrams can involve also other or- 
dered phases and states in addition to those obtained 
for Hamiltonian Q |40fH7] . even in the ground state. 
We can suppose that small but finite single electron hop- 
ping tij will not qualitatively alter the phase diagrams, 
at least for the case ksT > J^j tfj/U- The main effect of 
tij (for \U\ ^ tij , U < 0) is Si renormalization of the pair 
hopping term lij — ^ lij +t1j/\U\ and an introduction of 
an effective intersite density-density repulsion ~ tfj/\U\. 
For U <{) and / < the charge density wave state can 
also occur [251 SH miS] . For ?7 > and both signs of 
/ the tij term generates antiferromagnetic correlations 
(in particular for n = 1) [TJ I49| competing with super- 
conducting ones and external field B, and its effects can 
essentially modify the phase diagrams and the properties 
of normal state. In such a case it is necessary to consider 
also various magnetic orderings |4 H |42 | [50]. Moreover, 
several phase separation states involving superconduct- 
ing, charge and (or) magnetic orderings could also be 
stable for n 7^ 1. 

Although our model is (in several aspects) oversimpli- 
fied, it can be useful in qualitative analysis of experimen- 
tal data for real narrow-band materials with very short 
coherence length (exemplary systems mentioned below) 
and fermions on optical lattices [TBHin]- In particular, 
our results predict the existence of the electron phase sep- 
aration (SS-NO) and describe its possible evolution and 
phase transitions with increasing T and a change of n (or 

in the presence of external magnetic field. Notice that 
the temperature dependence of the upper critical field 
in unconventional superconductors has a positive curva- 
ture in coincidence with results of figure |4] Obviously 
such a PS state is different from the Abrikosov-Shubnikov 
mixed-state in type-II superconductors, e. g. no magnetic 
flux quantization, no vortex lattice, etc. Such „upper 
critical fleld" (at T = 0) is independent of n and depends 
only on [///q (it decreases with increasing U /Iq), whereas 
„lower critical fleld" decreases with increasing n for flxed 
U/Io- 

It is well known that there exist two distinct ways to 
induce pair-breaking in type-II superconductors by an 
applied magnetic fleld: orbital and spin-paramagnetic ef- 
fects. The former is related to an emergence of Abrikosov 
vortex lines. The spin-paramagnetic pair-breaking ef- 
fect comes from the Zeeman splitting of spin singlet 
Cooper pairs. The actual Hc2 of real materials is gen- 
erally influenced by the both these effects. The rela- 
tive importance of the orbital and spin-paramagnetic ef- 
fects can be described by the Maki parameter a [51J. 
For a ^ I (conventional superconductors) orbital effect 
dominates, whereas in materials with a heavy electron 
mass (narrow bands) or multiple small Fermi pockets 
a > 1 and paramagnetic effect becomes crucial. Some 
exemplary materials with a > 1 are, among others, 
CeCoIns, k-(ET)2Cu(NCS)2, a-(ET)2NH4Hg(SCN)4, k- 
(ET)2Cu[N(CN)2]Br and A-(BETS)2GaCl4. Notice also 
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that for interacting fermions on the non-rotating optical 
lattices only the paramagnetic effect can occur. 

In the present model the magnetic field only acts on the 
spin through the Zeeman term. An interesting problem 
of the orbital contribution through the pair hopping is 
left for future study. 

The phase separation instability is specific to the short- 
range nature of the model. The (unscreened) long- 
range Coulomb interactions prevent the large-scale PS 
of charged particles |48] and only a frustrated PS can oc- 
cur (mesoscale, nanoscale) with the formation of various 
possible textures [52 54 . 

The PS states involving superconductivity are shown 
experimentally in several systems. For example, or- 
ganic compounds exhibit the superconductor-insulator 
phase separations [6-8 , whereas mesoscopic phase sepa- 
ration has been observed in the family of iron-pnictides 
[T3j. Moreover, for special cases of La2Cu04+5 
and La2-a;Srj:Cu04+5, muon and superconducting quan- 
tum interference measurements suggest existence of fully 
phase separated regions [9UTT| . Finally, among the mate- 
rials for which the local electron pairing has been either 
established or suggested the best candidates to exhibit 
the phase separation phenomena are doped barium bis- 
mutates |14| [TSl |48] . Recent experiments on ultracold 
imbalanced Fermi gas trapped in external harmonic po- 



tential serve as an alternative way to study the pure Zee- 
man effect on Fermi superfluidity. Several evidences of 
phase separations in such systems to state containing a 
paired center core and unpaired atoms outside this core 
have been reported |l8ti20j . 

It is of interest to analyze the impact of density-density 
|37l [55] and magnetic |!55J |57] interactions on the phase 
diagrams of model Q. Some results concerning the in- 
terplay of these interactions with the pair hopping term 
for B = have been presented in (4j |50] . 
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